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ABSTRACT 

We extract the resonant orbits from an N- body bar tha t is a good representation of the Milky 
Way, using the method recently introduced by Molloy et ah (2015). By decomposing the bar into its 
constituent orbit families, we show that they are intimately connected to the boxy-peanut shape of the 
density. We highlight the imprint due solely to resonant orbits on the kinematic landscape towards 
the Galactic centre. The resonant orbits are shown to have distinct kinematic features and inay b e 
used to explain the cold velocity peak seen in the APOGEE commissioning data (jNidever et al.|2012|). 
We show that high velocity peaks are a natural consequence of the motions of stars in the 2:1 orbit 
family and that stars on other higher order resonances can contribute to the peaks. The locations of 
the peaks vary with bar angle and, with the tacit assumption that the observed peaks are due to the 
2:1 family, we find that the locations of the high velocity peaks correspond to bar angles in the range 
10° < ^bar < 25°. However, some important questions about the nature of the peaks remain, such 
as their apparent absence in other surveys of the Bulge and the deviations from symmetry between 
equivalent fields in the north and south. We show that the absence of a peak in surveys at higher 
latitudes is likely due to the combination of a less prominent peak and a lower number density of bar 
supporting orbits at these latitudes. 

Subject headings: Galaxy: kinematics and dynamics — Galaxy: evolution 


1. INTRODUCTION 

It is now widely accepted that the Milky Way (MW) hosts 
a bar. Many methods have been used to m ap out the struc¬ 
ture of the bar, such a s IR photometr y ([Blitz Spergel 
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Mess et al. |2U12|). bigniricant observational ettort has been 
expended in mapping the spatia l density of the bar using the 
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OGLE-III and VVV data (e.g., |Wegg fc Gerhard|[20T3 


Cao 


et al. | 2pi3[l. More recentl y, using a large sample ot red clump 


giants, I Wegg et al. (2015) derive a bar half-length of between 
~4.5 and 5 kpc. However, uncertainties about the nature 
of the interface between the bar and the disk or spiral arms 


can strongly in fluence these estimates (Martinez-Valpuesta <V 
Gerhard [2011 ). Uncertainty also surrounds some other tun- 
damental parameters of the bar, such as the viewing angle 
^bar and the pattern speed of its rotation Ubar- The litera¬ 
ture reports values for Oh^ir in the range ^ 20° < ^bar £ 45 


e.g., Stanek et al.||1997 Benjamin et al.|[2005 ) and for Ubar 


in the range ~ 25 < Qh, 

.p < 50 km s ^ kpc ^ (e.g.. 

Antoja 

et al.|2014 jPortail et al.j 

2015). 



Galactic bulge, kinematic substructure has rarely been ob¬ 
served. This is a pity, as such substructure may betray evi¬ 
dence of the processes that formed and shaped the bar. The 
Bulge RAdial Velocity Assay (BRAVA; Rich et al. 120071) and 


GIRAFFE Inner Bulge Survey (GIBS; |Z;Occali et al. 
both observed ~ 10,000 giants over a large region ot the 
but revealed no sig nature of cold high velocity peaks. The 


p014) 

bulge. 


ARGOS survey (see Ness &: Freeman 2012 Ness et aT||2013 ) 


has also yet to reveal evidence tor streams m tHe bulge, al- 
though their velocity distributions, cut according to metallic- 
ity, hint at the wealth of information contained in the kine¬ 
matic data. 

Recentl y, however, a cold hi gh velocity peak has been ob¬ 
served by Nidever et al. (2012) in the Apache Point Observa¬ 
tory Galactic Hvolution Experiment (APOGEE) commission¬ 
ing data. For certain fields towards the Galactic Bulge, they 
find bimodal velocity distributions and identify cold {a ~ 20 
km s“^) secondary peaks in the distribution of line of sight 
velociti es at uios. ^ 200 km s~L The independent observa¬ 


tions of Babusiaux et al. (2014) also hint at the presence of 
a high velocity peak, this time with red clump stars. The 
origin, and even the exis tence, of this feature has been th e 
subject of recent debate ( Li et al.]|2014 Zoccali et aL||2014 ). 
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Figure 1. The normalised surface density (contours at 0.1 dex) in the x-y (top) and x-z (bottom) planes for all (left) and 2:1 resonant 
orbits (right). 


Li et al.|(2014) found the absence of a statistically significant 
cold high velocity peak in two N-body barred models. They 
also showed that it is possible for a spurious high velocity peak 
to appear if there are only a limited number of stars observed. 
Here, we look at the matter anew, using novel algorithms to 
extract nearly periodic orbits from bar simulations. 

2. ORBITAL COMPONENTS IN THE BAR 


In Molloy et al. (2015 hereafter M15), we introduced a 
method to identity resonant orbits in N-body simulations, 
and used it to provide a possible explanation for the bi- 
modal velocity distributions observed towards the Galactic 
anti-Centre. The method is an alternative to identifying pe¬ 
riodic orbits using frequency analysis and is well suited to 
analysis of evolving N-body models with a changing poten¬ 
tial (see M15 for details). It is now applied to the inner parts 
of a barred model of the Milky Way. Here, a cold but thick¬ 
ened disc self-consistently develops a bar, which undergoes a 
buckling instability to forn i a Bulge that is a good match to 
that of the Milky Way (see Shen et al. pMoj [!T¥'Shen|2012| 
Qin et al. 20I5|). Resonant orbits can be characterized by 
the tact that they close and return to a previously occupied 
location in phase space in some rotating frame. The method 
in MI5 proceeds by recalculating the N-body orbits in many 
different rotating frames. We define a metric Dps to measure 
the distance each particle travels from some arbitrarily chosen 
starting point in the rotating frame. If an orbit almost closes. 
Dps should, at some point along its trajectory, be nearly zero. 
By defining a cut-off, we can extract a sample of the nearly 
closed orbits from the simulation. Resonant orbits can librate 
about each family’s parent orbit, so by defining tighter and 
tighter cuts on Dps, we can extract cleaner and cleaner sam¬ 
ples of resonant orbits. The choice of cut on Dps is really set 


by the problem in hand. 

We make some minor modifications to the algorithm de¬ 
scribed in MI5. In bars, there are always some chaotic or¬ 
bits, especially near corotation. Indeed, chaotic orbits may 
make significant contributions to the structure and dynamics 
in barred systems. It has been shown that the fra ction of 
chaotic orbits is sensitive to the bar strength and size (Manos 
& Athanassoula|20II| |Manos & Machado|20I4|). Chaotic or¬ 
bits may return arbitrarily close to their chosen starting point 
over long timescales (the Poincare Recurrence Theorem), and 
can therefore be mistaken as periodic. Previously, we mea¬ 
sured the phase space distance from a single point as the 
orbit proceeded on its trajectory over ~ I Gyr. Here, we 
define a time frame, unique to each particle, over which we 
apply the phase space distance method. For each particle, we 
measure the duration it takes to complete eight radial oscil¬ 
lations. This ensures that we have a long enough trajectory 
to extract high-order periodic orbits while excluding chaotic 
orbits that rapidly explore their phase space volume. For or¬ 
bits very close to the centre, this time frame may be sampled 
by very few points so we interpolate the trajectory. 

Also, instead of measuring the phase space distance from 
a single point, we find the time at which a particle reaches 
its first apocentre, to- As we scan different rotating frames 
(between 37 < Qp < 40 km/s/kpc), we extract sections of 
the trajectories that lie in the range (f)'{to) to (f)'{to) + 7r/4, 
where <j)' is the azimuthal angle in the rotating frame. This 
gives us, say, n similar sections of the trajectory. If n is 
less than three, then we increase the duration over which we 
apply the method. The n — 1 section is the reference section 
which we compare to following test sections. The phase space 
distance Dps is measured between successive points on the 
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q = (Q-QJ/k 


Figure 2. The distribution of q = {Q. — Q^c )/for the sample of 
periodic orbits with Rg < 2.5 kpc. We select 2:1 orbits as those 
with 0.48 < q < 0.52, 3:1 orbits as those with 0.33 < q < 0.35 and 
5:2 orbits as those with 0.38 < q < 0.42. 


reference section and the test sections. We take the average 
Dps between the orbit sections as a measure of how “closed” 
the orbit is - the lower the value, the more closed the orbit. 
Of the n — 1 Dps values, we take the minimum and in the 
following we adopt a cut on the phase space distance of Dps < 
0.06 (Over the course of an orbit Dps varies between 0 < 
Dps < a / 2 , see M15). 

By extracting a sample of resonant orbits in the central 
regions, we can deduce the contribution they make to the 
structure of the bar. Figure shows the normalized sur¬ 
face density of the inner parts of the disc. In the left pan¬ 
els, we show the x-y (top) and x-z (bottom) surface den¬ 
sity for all of the particles in the simulation. The bar ex¬ 
tends to ~4 kpc or Rck/cl ~ 1.125 and has an axis ratio of 
h/a ^ 0.5. This is s lightly less extended th an the value of 
0.35 found by OGLE (|Rattenbury et al.|2007|) but agrees well 
with t he structure derived rece ntly trom a large sample RR 
Lyrae (|Pietrukowicz et al.|2014|) . The x-z surface density ex¬ 
hibits a strong boxy-peanut shape characteristic of buckled 
bars. We extract the resonant 2:1 orbits by estimating the 
azimuthal (Q) and epicyclic (k) frequencies for our sample of 
closed orbits. For each orbit, we calculate g = (Q — Qp)//^, 
where Qp is the pattern speed of the frame in which the or¬ 
bit closes, or reaches its lowest Dps. We then extract the 2:1 
orbits as those having 0.48 < q < 0.52. The surface densities 
are shown in the right panels of Figure It is clear that the 
2:1 orbits generate the backbone of the bar. The buckling 
instability has a noticeable effect on this family of orbits, in¬ 
ducing a large vertical velocity dispersion for stars at the end 
of the bar. The contribution of the 2:1 family to the boxy 
and peanut shape is unmistakable. 

Of course, other orbit families are present in the simula¬ 
tion. Figure!^ shows the distribution of ^-values for particles 
with Rg < T5 kpc. The most prominent family are the 2:1 
orbits but there is also a significant contribution from the 3:1 
{q ~ 0.33) and 5:2 {q ^ 0.4) orbit families. The relative frac¬ 
tions for these families is approximately 2:1:1 (where we have 
taken the 3:1 orbits as those with 0.31 < q < 0.35 and the 
5:2 orbits as those with 0.38 < q < 0.42). The relative frac¬ 
tions between the orbit families remains largely unchanged 
using different cuts on Dps- However, the normalisation of 
the phase space distance (see Ml5) affects the elongated 2:1 
orbit family differently compared to the more circular 3:1 and 
5:2 orbits. Consequently, we refrain from making any strong 
conclusions based on the relative fractions of the extracted 
periodic orbits. In any case, we choose a cut on the phase 
space distance that allows us to sample the phase space of 
each family sufficiently. So, while the relative fractions of 
each family is uncertain, the shapes of the distribution func¬ 
tions are robust. 


3. VELOCITY DISTRIBUTIONS 

The photometry and star counts for this model have already 
been shown to be a good match to observations. Indeed, the 


Table 1 

Line-of-sight velocity distributions for 
each of the APOGEE commissioning 
fields (truncated). Table 1 is published in 
its entirety in the electronic edition of 
ApJ, a portion is shown here for guidance 
regarding its form and content. 


(l,b) 

^bar 

All Peak^ 

2:1 Peak^^ 

n 


km s“^ 

km 


10^ 

60.2 

192.9 

(10.0,-2.0) 

15° 

20° 

55.2 

56.3 

192.4 

175.6 


25° 

51.9 

180.9 


10“ 

70.1 

197.3 

(10.0,2.0) 

15° 

20° 

64.7 

62.5 

193.9 

197.7 


25° 

56.3 

183.9 


^ ffis to all particles in the field. 
^ Maximum peak of 3-Gaussian fit. 


simulation was tailored to match the kine matics towards th e 
Galactic Bulge as seen by the BRAVA data (|Shen et al.|2010|). 
The pattern speed of the simulation bar stays roughly con- 
stant at ~40 km s“^ kpc“^ (see Ml5). This is in the middle of 
the range of values reported in the literature and is c onsistent 


with the m ost recent estimates from gas dynamics (Sormani 
et al.||2015| p" A suite of simulations spanning the parameter 


space ot baf size, strength, viewing angle and pattern speed 
is required to fully investigate their effect on the r'los distri¬ 
butions. 

Since we can deconstruct the bar into its different orbital 
families, we can now characterize the contribution of each 
family to the velocity distributions. The APOGEE commis¬ 
sioning data re vealed cold high vel ocity peaks for a number 
of Bulge fields. Binney et al. (1991) were the first to link the 
motions of gas towards tfie Bulge to orbits in (planar) ellip¬ 
tical potentials. Below, we will show that these peaks arise 
naturally as a result of the motions of resonant bar orbits, in 
particular, the 2:1 orbital family. Note that if we include all 
simulatio n particles in o ur v\os distributions we recover the 
results of Li et al. (2014), in which no cold peaks are found. 

To generate mock v\os distributions, we first fix some funda¬ 
mental parameters. We choose the Solar radius as Ro = 8.5 
kpc and the circular velocity at Rq as Vc = 220 km s“^. 
Varying these between reasonable values has only minor ef¬ 
fects on the distributions. We assume an angle between the 
long axis of the bar and the Solar-Galactic Gentre (GG) line 
of ^bar = 15° (we later justify the choice, where we explore 
a range of bar angles). We also limit the distances of the 
particles to between 3 kpc and 9 kpc and, in order to in¬ 
crease the numbers in the samples, we include particles from 
equivalent positions on either side of the disc and increase 
the diameter of the field by a factor of two compared to the 
APOGEE fields. As a further measure to increase the number 
of particles, we also average over 10 timesteps, making sure 
to take into account the (small) change in bar angle between 
timesteps. 

The top panel of Figure shows t he kinematic data on 
the APOGEE commissioning fields of Nidever et al. (2012 
together with their two Gaussian decomposition, in the bo 
tom plot of Figure we show the velocity distributions for 
our sample of 2:1 orbits. To avoid forcing fits to binned 
data, we instead opt for a more general approach. We pop¬ 
ulate Gaussian mixture models (GMMs) using a genetic al¬ 
gorithm that converges on probability distribution functions 
that could have produced the datc[^ This only requires one 
to input the range of parameter space to explore. In order 


}: 


^Aumer & Sclionrich| (|2015| have recently shown that the vios 
distributions are only weakly atiected by the pattern speed, at least 
between 25 and 30 km s“^ kpc“^. 

^We do this using the freely available SOLDER routines: 
http:/ / WWW. ast.cam.ac.uk/~vasily/solber/ 
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Figure 3. Top: For convenience, we show that distributions from |Nidever et al.| (|2012|). The blue and red curves show their 
decomposition of the velocity histogram into two Gaussians. Bottom: The i;ios distributions tor z:l resonant orbits in each field, assuming 
a bar angle of 15°. A genetic algorithm has been used to populate a 3-Ganssian m i xture model for the data. The underlying Gaussians 
are shown as blue curves and the locations of the peaks reported in |Nidever et al.| (|2012|) as the vertical dashed lines. The positions of 
the peaks show the locations at which 2:1 resonant orbits make the largest contribution to the observed histograms, which of course also 
include contributions from other orbital families. 
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Figure 4. A schematic (and simplified) representation of a 2:1 bar orbit. The red arrows represent the motion of the particle in the frame 
of the bar (the black dashed box). For different fields, represented by the blue dashed lines, it is clear that this type of orbit gives different 
contributions to the line of sight velocities. For / > 0°, we expect the star to be moving away, whilst the reverse is true for fields with 
/ < 0°. Inset, top right: A sample orbit from a high velocity peak at (l^b) = (10,0), the red triangles indicate the fields of view from 
equivalent sides of the disc. The blue dots represent the preceding 100 timesteps of the orbit and the large blue cross, inside the field of 
view, indicates the position of the particle when it has a high v\os (listed). Inset, bottom right: The density of the 2:1 resonant orbits 
in the field at (/, b) = (10, 0) as a function of v\os and distance. The highest velocity particles lie at a distance of ~ 8.5 kpc corresponding 
to pericentric passage. The lowest velocity particles lie at a distance of ~ 6 kpc corresponding to apocentric passage. 


to compare the fits across each field, we force the distribu¬ 
tions to be fit with three Gaussians. Initially, we fit each of 
the field’s distributions with one to five Gaussians. We then 
performed likelihood ratio tests to see when adding an extra 
Gaussian component made no significant improvement. Most 
of the fields preferred either two (43.75%) or three (52.5%) 
Gaussians, while only a small proportion preferred one Gaus¬ 
sian (3.75%). Generally, the distributions are split into neg¬ 
ative and positive velocity components, with an intermediate 
component in some fields. We interpret the negative velocity 
component as being due to particles on the near side of the 
bar, streaming towards apocentre. The high velocity compo¬ 
nent are then the particles streaming towards pericentre on 
the far side of the bar, while the intermediate component rep¬ 
resents the particles that are slowing down as they approach 
apocentre, those at apocentre (with almost zero line-of-sight 
velocity) and those leaving apocentre, speeding up as they 
head towards pericentre. The shape, and number of compo¬ 
nents in the distribution is a non-trivial function of the field 
being observed, the chosen bar angle and the range of dis¬ 
tances being sampled. For a selection of bar angles, we list 
the values of the peaks of the distributions for each field in 
Table [2 In the case of the three Gaussian fits to the 2:1 
distributions, we list the highest valued peaks. 

A crude representation of a 2:1 bar orbit is shown in Figure 
1^ The orbit, shown as the red arrows, reaches its apocentre 
at each end of the bar (on the x-axis), while the pericen- 
tres are located in the directions perpendicular to the bar 
(the black dashed box). The black dot represents the posi¬ 
tion of the Sun, so that the bar is rotating in a clockwise 
direction and the Galactic Gentre-Solar position line is at an 


angle close to 20° with respect to the bar. The dashed blue 
lines indicate the fields of view in Galactic coordinates. With 
(/, b) = (10°, 0°), we can see how this type of orbit contributes 
to a high velocity component in this direction. For this line- 
of-sight, the measured uios captures most of the components 
of Galactocentric ur and and results in the high velocity 
peaks shown in Figure!^ In the Galactocentric frame, these 
stars have ur < 0 km s~ and ^ 0 km s“^ which, for this 
particular direction in relation to the Solar position, make 
a significant contribution to uios “ in fact, almost all of the 
particle’s velocity is coincident with the line of sight. 

As a specific example, we take a sample orbit in the di¬ 
rection of (/,5) = (10°, 0°) from the high velocity component 
('^^los >150 km s“^) in the velocity distributions assuming a 
bar angln of 20°. We plot the orbit in the upper right panel of 
Figure]^ and indicate with the red lines the fields of view from 
which our sample is derived. The orbits are plotted over a 
period of ~1 Gyr with the final 100 timesteps indicated with 
blue crosses and the final timestep shown as the large blue 
cross lying inside the field of view. The line of sight velocity, 
assuming a Solar position of {x,y) = (8.5,0), is also listed. 

A population of near-circular orbits would make a similar, 
but smaller, contribution to the high velocity peak in uios 
since ur 0 km . This direction also captures stars that 
are just reaching their maximum radius and so contribute to 
the negative velocity component in the distributions of uios- 
The peaks of this component are at a lower velocity compared 
to the simple disc rotation model. This shows that resonant 
orbits in this direction imprint both low and high velocity 
kinematic signatures on the line of sight velocity distribu¬ 
tions. A field centered on (/,5) = (0°,0°) passes through the 
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Figure 5. The ^’los distributions for the 3:1 (top two rows) and 5:2 (middle two rows) orbit families in the APOGEE commissioning fields. 
For the 3:1 and 5:2 families we average over 35 & 50 timesteps (~0.34 & 0.48 Gyr) respectively to generate the dis tributions, which ar e 
then modeled with GMMs. The location of the observed peak for each field is indicated by the vertical dashed line ( [Nidever et al. |2012| >. 
The distributions due the 2:1, 3:1 and 5:2 families combined are shown in the bottom two rows (where we have averaged over iO tirnesteps 
or ~100 Myr). The relative fractions of each family are listed inset (2:1-red, 3:1-blue, 5:2-green). We overplot our simple model which uses 
relative fractions of 2:1:1 (for the 2:1, 3:1 and 5:2 families respectively) as the black curve. The components, using these relative fractions, 
are shown as the coloured curves where we use the 3 Gaussian fit to the 2:1 family from Figure and the fits to the 3:1 and 5:2 families 
from the top four panels. 
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Figure 6. Left Panels: For the APOGEE commissioning fields shown in |Nidever et al. 
characterize the vi^s distributions assuming bar angles in the range 5^ 


we use Gaussian mixture models to 
< ^bar <' 5U^ (shown here is ill ° < ^bar < 25°). The black stars/dots 
represent the general/high velocity peaks seen in the APOGEE data (note that some longitudes have two points since APOGEE observed 
fields symmetric about the plane). The blue dots are the highest velocity peaks from our Gaussian mixture models for the 2:1 orbit velocity 
distributions. The red stars are the peaks from fifs to all particles in the fields. Right Panel: The difference in the values of the 
simulation and observed peaks for fields with I < 10° as a function of ^bar, to which a linear fit is made (black line). The intersection of 
the fit to the 0 km/s difference suggests a bar angle of ^bar = 17.7°. 


whole structure of the bar and therefore catches stars with a 
negative uios on the near side and with positive v\os on the 
far side. At negative /, the bar stars are approaching the 
Solar position and so imprint a high negative velocity compo¬ 
nent, mirroring the corresponding fields at positive 1. Having 
a non-zero bar angle influences the differing shapes of the 
distributions between positive and negative longitude fields. 

The lower right panel of Figure]^ shows the density of par¬ 
ticles as a function of uios and distance. For the line of sight 
(/,5) = (10°, 0°), particles that are reaching apocentre con¬ 
tribute to the peak at ~ 0 km s“^ and are at smaller dis¬ 
tances (~6 kpc). As the distance is increased, the particle’s 
velocity increases since both and ur are increasing and 
also because the velocity vector is coincident with the line of 
sight. Of course this is only true for the 2:1 orbits, the wide 
range of other orbit families in static and evolvi ng potentials 


will have more complicated morpholo gies (e.g., Pfenniger & 
Friedli||1991| |Manos & Machado||2014| . The highest velocity 
particles are approaclimg their pericentre occurring on the 
short axis of the bar at a distance of ~8 kpc, very close to 
Rq. This picture is consistent with our interpretation of the 
high velocity peaks in Bulge fields being due to the motions 
of resonant 2:1 bar orbits. 

The other orbit families are also likely to have an influence 
on the uios distributions. We see from Figure that the 3:1 
and 5:1 families make up a significant portion or our periodic 
orbits. If we plot the uios distributions for these orbits it is 
clear that they are rich in structure (Figure [^. To account 
for their lower numbers, we average over 35 arm 50 timesteps 
(336 Myr and 480 Myr) for the 3:1 and 5:2 families respec¬ 
tively. In this way, each of the chosen particles moves around 
their respective orbit pattern roughly twice. For fields with 
I > 10° the 3:1 (top two rows) and 5:2 (middle two rows) fam¬ 
ilies also generate strong high velocity peaks (as before, we 
use a likelihood calculation to choose the number of Gaus¬ 
sian components). Just as for the 2:1 families, the peaks 
in these fiel ds he at a lower veloc ity compared to the peaks 
identified by|Nidever et al.|(2012) (vertical dashed lines). At 
lower longitrrdes (/ < t)°) these tamilies also generate peaks 
coincident with the peaks identified in the data (e.g., 3:1 and 
5:2 at (/,5) = (5.7°,—2.0°)). However, there are fields in 
which the peaks fail to match the data (e.g., 3:1 and 5:2 at 
(/, b) = (4.0°, 0.0°)). For certain bar angles ^bar the 2:1 family 
can simultaneously, across a number of fields, generate peaks 
that are a good match to the data. 

The relative fractions of each orbit family is an important 
factor in shaping the uios distributions (i.e., the locations of 


cold peaks). It is likely the case that the relative fraction in 
each field differs from the “global” fraction for the bulge as 
a whole. This is because each family has a distinct spatial 
density distribution. For example, the 2:1 family stays close 
to the long axis of the bar while the other families generally 
travel farther along the short axis. This, together with the 
uncertain distances probed in the APOGEE fields, makes it 
very difficult to derive the relative fractions from the data. 
We show the uios distributions imprinted by the 2:1, 3:1 and 
5:2 families combined in the bottom two rows of Figure 
The combined 2:1, 3:1 and 5:2 uios distributions average over 
10 timesteps. As a fiducial model we apply simple relative 
fractions to the fits for each family (2:1:1 for the 2:1, 3:1 
and 5:2 families respectively; for the 2:1 family we use the 
3-Gaussian fits from the previous section) and overplot the 
simple model on the uios distributions. The actual fractions 
from the model, listed in each panel, vary only slightly from 
field to field but for this exercise we have applied a uniform 
distance cut along each field. The distances probed in the 
data are likely to vary significantly with Galactic I and b. 
Given the small change in fractions between fields, it’s not 
surprising that the simple model matches the distributions 
quite well. With this model the high velocity peak is largely 
determined by the 2:1 family, especially for the fields with 
5 = 0°. At larger latitudes, the other families appear to make 
more significant contributions to the high velocity component, 
but never stronger than the 2:1 family. 

With the tacit assumption that the locations of the peaks 
are set by the 2:1 family (i.e., that the relative fractions of 
the most populated families is close to 2:1:1, as above), can 
we use this insight to constrain the viewing angle of the bar, 
using the information on how the velocity peaks vary with 
Galactic position? In Figure (left panels), the bl ack stars 
and dots rep resent the locations of the peaks found by |Nidever| 
et al.| (2012). Over-plotted as blue dots and red stars are the 


positions ot the peaks from our models. Specifically, the blue 
dots are derived from the Gaussian mixture models for the 
2:1 resonant orbits, whilst the red stars are extracted from 
fits to ah the particles in the field. It is clear that, although 
no single choice of viewing angle reproduces all the data, the 
trends in the velocity peaks are well-reproduced for bar angles 
- 15°. 

The fields with the largest deviations are the ones with 
the highest longitude, for which our simple picture probably 
breaks down. Although the peaks in these fields are quite 
pronounced, they systematically he at lower values compared 
to the data. This could be due to one of two possible scenar- 
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Figure 7. Top Left: The projected surface density of the 2:1 resonant orbits assuming a bar angle of ^bar = 15° (averaged over 0.48 
Gyr). Horizontal lines indicate stripes in longitude from which the following histograms are drawn. (A) uiog distribution for the strip with 
— 10° < I < 10° and 6 = 0°. The black histogram represents a strip of width Ah = 1° (thick) and the red, a strip of width Ah = 0.5° 
(thin). (B) As for panel A, but at 6 = —2°. (C) As for panel A, but at 6 = —4°. We also make fils to the distributions, the mean and 
dispersion of which are listed for the thick and thin strips respectively. 


ios. Firstly, being at high longitudes, particles in these fields 
feel a significant cumulative effect of the shallower potential. 
The high velocity particles in fields with I > 10° are at peri- 
centre between 1.5 and 2.0 kpc along the short axis of the bar. 
We expect the potential to be somewhat shallower since the 
pure disc simulation is absent of a gaseous component and 
live halo that may relax into a more concentrated configura¬ 
tion after the formation of the bar. However, the simulation 
has been shown to be in good agreement with the kinematics 
observed by BRAVA (|Shen et alJl^lOj, even as far out as 

To 


I = 10°. A good match is made to the mean velocities and 


velocity dispersions, so that the comparison is made to the 
data through the whole line of sight. As we’ve shown above 
(Figure® bottom right), the high velocity peaks correspond 
to a limned range in distance. That the deviation from the 
observed peaks increases with I is another indication that a 
somewhat shallow potential is the cause. 

Another possible scenario is that the peaks are in fact 
caused by another family of resonant orbits. We have checked 
the other major families in the bar, the 3:1 and 5:2 families. 
They do generate strong peaks in these fields but, as with the 
2:1 orbits, the high velocity peaks are at systematically lower 
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values (combinations of the 2:1, 3:1 and 5:2 orbits also result 
in peaks with low values, see Figure!^. However, other higher 
order resonan ces may also be impo ^ant in these regions. As 
mentioned in Nidever et aL| (2012|, the high velocity peaks 
are unlikely to be due to tidal streams. Although the Sagit¬ 
tarius stream lies close by on the plane of the sky, the high 
velocity stars show no preference for magnitude or metallicity 
and, in any case, t he stream stars are not expected to appear 
in large numbers (Law & Majewski|2010 ). 

Another way of synthesizing this mtormation is presented 
in the right panel of Figure [m which shows how the differ¬ 
ence between the observed arm simulation high velocity peak 
varies with assumed viewing angle of the bar. The observa¬ 
tional data for each field are represented by the horizontal 
red dot-dashed line (i.e., the zero-difference line), and sug¬ 
gest viewing angles between 10° and 25°. Some scatter is 
expected, as the N-body model does not exactly reproduce 
the three-dimensional density of the inner Galaxy and the 
kinematical properties are subject to numerical shot noise. 

4. DISCUSSION 

We have shown that the high velocity peaks seen in the 
APOGEE commissioning data may be explained by the pres¬ 
ence of a large family of 2:1 resonant orbits in the Galacti c 
bar, as was tentatively suggested by Nidever et al. (2012|). 
These orbits are elongated along the bar's major axis, and 
must support the distended shape of the bar t o provide its 
backb one. Indeed, it has already been shown by|Binney et al.| 
(1991) that the motion of gas towards the Bulge follows nat¬ 
urally from orbits in (planar) ellipsoidal potentials. When 
viewed at bar angles in the range 10° < ^bar < 25°, the 2:1 
orbits naturally give rise to secondary peaks in the line of 
sight velocity distributions at uios ~ 200 km s“^. We have 
provided a pictorial explanation of this phenomenon. 

Our interpretation is open to the objection that the method 
is not fully self-consistent. We have shown that the popu¬ 
lation of 2:1 orbits can generate t h e kin ematic features to 
explain the data of Nidever et al. (2012), but we have al¬ 
lowed the normalization ot the density m these orbits to vary 
independently of the N-body model from which they were ex¬ 
tracted. This though is unlikely to be a serious concern, as 
the range of self-consistent equilibria for bars is wide and so¬ 
lutions will exist using different relative populations for the 
orbital families that comprise the bar. This has been demon¬ 
strated expli citly for the related problem of static triaxial el¬ 
lipsoids (e.g.,|Statler|1987||Hunter|1995|. Eor bars, the ability 
to exchange orbits without changing the density is indicated 
by both the non-unique decom position into the clas sical and 
regular components found by |Hafner et al. (2000| and the 
range of made-to-measure m odels reproducmgf the d ensity in- 
ferred from the VVV survey (|Wegg & Gerhard|2013| found by 
4s we tiave mentioned above, t 
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the other 


resonant orbit tamilies also produce rich structure in their ve¬ 
locity distributions. If we combine the 2:1, 3:1 and 5:2 orbits 
in our distributions, the high velocity peaks (which are some¬ 
what more obscured) are dominated by the 2:1 orbits and so 
follow the trends outlined above to suggest ^bar ~ 15°. An¬ 
other important consideration is the range of distances probed 
by observations. Given the non-uniform distribution of dust 
towards the Galactic Bulge, it is certain that different fields 
are reaching different distances. Indeed the range of distances 
reached may vary significantly through just one field. The 
lower right panel of Eigurej^ shows that the highest velocity 
2:1 stars are placed at a distance close to Rq. There is also 
likely to be some variation in the relative fractions of the dif¬ 
ferent orbit families as one varies the distance probed. It is 
clear then that distance has an important role in shaping the 
velocity distributions. 

4.1. Puzzles in the Data 

Some issues remain about the high velocity peaks. The first 
is why symmetry is not seen between positive and negative 
latitude fields. Peaks are observed in the field at (4.3,-4.3) 
but not in the field at (4.3,4.3), the same goes for the fields 


at (5.7,-2) and (5.7,2). The obvious explanation for such a 
difference is that the distances being probed differs between 
north and south. It is kn own that extinction in the north is 
greater than in the south (Gonzalez et al.|2012), which offers 
a possible explanation for tne disparity. However, according 
to APOGEE estimates, although the average extinction is 
indeed greater for the field at (4.3,4.3) than for the field at 
(4.3,-4.3), the opposite is the case for the fields at (5.7,2) and 
(5.7,-2). The extinction data is at its most uncertain close to 
the plane of the disc, so it is unclear if this is causing the dif¬ 
ference between symmetric fields. Another point to consider 
is that in the plane, where extinction is highest, strong peaks 
are seen. As we show below, a likely explanation here is that 
the number density of bar supporting orbits is higher in these 
regions. 

The second major concern is the differences between differ¬ 
ent Bulge surveys. The BRAVA, ARGOS and GIBS surveys 
don’t report the detection of cold components in their uios 
distributions. The presence of the peaks in the APOGEE 
data should be resolutely confirmed on analysis of the post¬ 
commissioning data. That the cold peaks are seen in one 
survey and not the others provides a possible clue as to the 
nature of bar supporting orbits. The clue lies in the dif¬ 
ferent observing strategies and selection functions for the 
data. The color selections are similar, with each survey using 
cuts on J — K color, but the magnitude ranges are not the 
same. These surveys each cover different footprints, with only 
APOGEE and GIBS observing a significant number of fields 
below \b\ = 4°. Most of our sample of resonant 2:1 orbits are 
confined close to the plane with \b\ < 5° (corresponding to 
~0.75 kpc at Rq). Eigure (top left) shows the projected 
surface density (on a log scale; contours at 0.1 dex) for the 
2:1 family, where we have averaged over 0.48 Gyr. In the 
following panels we show the uios distributions for stripes in 
longitude (—10° < I < 10°) centred on 5 = 0°,—2° V —4° 
(A, B V G respectively; black histograms indicate a strip of 
width Ab = 1° and red, a strip of width Ab = 0.5°). There 
is a stark difference between the distribution at 5 = 0° (A) 
and b = —4° (G). The strong peaks, clearly visible in the 
mid-plane, rapidly decrease in prominence as latitude is in¬ 
creased. At b = —4° there is only a small hint of a high 
velocity component at ~ 200 km/s. This, along with the low 
density of 2:1 orbits, makes the detection of a cold peak much 
less likely at higher latitudes. It seems reasonable then that 
the BRAVA and ARGOS surveys don’t report the detection 
of cold peaks. This effect is unlikely to be the only impor¬ 
tant difference between the surveys. As we have shown above 
(Eigure[^, the distances probed will also have a significant 
impact on the shape of the uios distributions. 

A more subtle difference is that the surveys observe dif¬ 
ferent types of stars (GIBS V ARGOS: mainly red clump; 
APOGEE V BRAVA: mainly M-giantsQ Since the targets 
that make up each survey differ in spatial (survey footprint, 
distance), temporal (ages) and chemical (color, magnitude) 
attributes, the question of why cold peaks are seen in one sur¬ 
vey and not others is certainly challenging (it should also be 
noted that the observing strategy may be selecting stars with 
a low velocity dispersion, thereby thinning the overall distri¬ 
butions to reveal the cold peak). A better question might be 
whether APOGEE has stumbled on a selection strategy that 
preferentially selects bar stars, and if so, can the strategy be 
shown to be consistent with chemical and dynamical models 
of the Galactic Bulge? 

We have so far not discussed in detail the effect of fore- 
/background or cha otic stars on the uios distributions. |Aumer 
& Schonrich (2015) suggest that the fore-/background stars 
play an important role in shaping the “main” distribution. 
The method outlined in M15 currently has no way of reliably 
identifying chaotic orbits but, in any case, they are unlikely to 
generate strong, long-lived features. Debattista et al. (2015) 


^It has been subsequently shown by|Aumer & Schonrich|(|2015| 
that the APOGEE selection function is indeed selecting young 
stars. 
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have also recently suggested that substructure, similar to that 
seen in the APOGEE data, can be generated by a rotating 
nuclear disc. Their kpc-scale nuclear disc provides an alter¬ 
native explanation for the observed kinematics. 

So far, we have considered t he cold peak to be a real and 
statistically significant feature. |Li et al.| (|2014| showed that 
peaks can be generated by under-sampling the wings of a 
broad r'los distributions. It is unlikely however that similar 
features wo uld be randomly drawn in a number of indepen¬ 
dent fields. Debattista et al. (2015) have also recently shown 
that such features are unlikely to arise due to Poisson noise. 

5. CONCLUSIONS 

We have set out a framework for understanding the sub- 
structure seen in th e APOGEE data. As was suggested by 
Nidever et al. (2012), we have shown that the motions of stars 
along periodic bar orbits naturally give rise to cold peaks in 
t^ios distributions. In particular, the 2:1 family are able to 
simultaneously generate peaks close to the observed values 
across a number of fields. This is a strong argument for the 
2:1 family being the main driver behind the location of the 
peaks. Although uncertain, the relative fraction for the differ¬ 
ent f amilies suggested by this s imulation (and previous works, 
e.g.,|Sparke & Sellwood||1987| points to the 2:1 family being 
the most populous m barred potentials. 

With this hypothesis in mind, we can constrain the view¬ 
ing angle of the bar with the locations the peaks generated by 
the 2:1 orbits to be 10° < ^bar ^ 25°. Although there remain 
uncertainties attached to this estimate, we are confident that 
large viewing angles can be discounted. The effect of funda¬ 
mental bar parameters on the locations of the peaks remains 
to be explored in any detail. Uncertainties associated with 
the stellar distances, the effects of extinction and the selection 
function of the APOGEE survey make a complete comparison 
with models difficult . We h ave also shown, in agreement with 
Aumer & Schonrich| ( 2015| , that the higher order resonances 
can make a sigmticant contribution to the r'los distributions, 
especially at high latitudes {\b\ > 2°). Glose to the plane, 
however, the 2:1 family appears to play the dominant role 
and sets the locations of any high velocity peaks. 

Einally, this work has demonstrated the power of the 
method introduced in Molloy et al. (2015) for extracting peri¬ 
odic orbits from N-body simulations m which the underlying 
gravitational potential is unsteady. By dissecting the bar into 
its constituent periodic families, we are able to study their 
characteristic kinematical and spatial signatures. Eorthcom- 
ing wide field spectroscopic surveys of the Galactic centre are 
likely to discover further streams, kinematic features and sub¬ 
structure. We are confident that the techniques of the paper 
will have an important role to play. 
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